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Abstract 



The Generalized Traveling Salesman Problem (GTSP) is an extension of the well-known Traveling Salesman Problem 
(TSP), where the node set is partitioned into clusters, and the objective is to find the shortest cycle visiting each cluster 
exactly once. In this paper, we present a new hybrid Ant Colony System (ACS) algorithm for the symmetric GTSP. 
The proposed algorithm is a modification of a simple ACS for the TSP improved by an efficient GTSP-specific local 
search procedure. Our extensive computational experiments show that the use of the local search procedure dramatically 
improves the performance of the ACS algorithm, making it one of the most successful GTSP metaheuristics to date. 
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1. Introduction 

The Generalized Traveling Salesman Problem (GTSP) 
is defined as follows. Let V = {1, 2, . . . , n} be a set of n 
nodes being partitioned into m non-empty subsets V = 
Ci U C 2 U . . . U C m called clusters. Let C(v) = Ci if v G Q. 
We are given a cost d uv of travelling between two nodes u 
and v for every u,v £V such that C(u) ^ C(v). Note that 
we consider only the symmetric case, i.e., d uv = d vu for any 
U, v G V, C(u) ^ C(v). Let T be an ordered set of nodes of 
size to such that C(Ti) ^ C(Tj) for j 6 {1,2,..., to}. 
We call such a set tour, and the weight of a tour T is 



w{T) 



i=l 



J Ti,Ti. 



(1) 



The objective of the GTSP is to find a tour T that mini- 
mizes w(T). 

It is sometimes convenient to consider the GTSP as a 
graph problem. Let G = (V, E) be a weighted undirected 
graph such that (u,v) E E for every u,v G V if C(u) ^ 
C(v). The weight of an edge (u, v) is d uv . The objective 
is to find a cycle in G such that it visits exactly one node 
in Ci for i = 1, 2, . . . , m and its weight is minimized. 

As a mixed integer program, the GTSP can be formu- 
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lated as follows: 
Minimize ^ 

(u,v)£E 

subject to 

^ x uv = ^2 x vu = y v for v eV, 

(u,v)£E (u,v)eE 



for i = 1, 2, . . . , to, 



z u — z v + (to — l)x uv < to — 2 for (u, v) G E, 

u>+ i,v i= l, 

x uv G {0, 1} for (u,v) e E, 

1 < z v < m — 1 for v G V \ {1}. 

The GTSP is an NP-hard problem. Indeed, if \d \ = 1 
for % — 1,2,..., to, the GTSP is reduced to the Traveling 
Salesman Problem (TSP). Hence, the TSP is a special case 
of the GTSP. Since the TSP is known to be NP-hard, the 
GTSP is also NP-hard. 

The GTSP has a lot of applications in warehouse or- 
der picking with multiple stock locations, sequencing com- 
puter files, postal routing, airport selection and routing for 
courier planes, and some others, see, e.g., 0] and references 
therein. 

Much attention was paid to the question of solving the 
GTSP. Several researchers proposed transformations of 
a GTSP instance into a TSP instance, see, e.g., [l[. At 
first glance, the idea of transforming a little-studied prob- 
lem into a well-known one seems to be promising. How- 
ever, this approach has a limited application. Indeed, 
such a transformation produces TSP instances where only 
the tours of some special structure correspond to feasible 



GTSP tours. In particular, such tours cannot include cer- 
tain edges. This is achieved by assigning large weights to 
such edges making the TSP instance unusual for the exact 
solvers. At the same time, solving the obtained TSP with 
a heuristic that does not guarantee any solution quality 
may produce a TSP tour corresponding to an infcasible 
GTSP tour. 

A more efficient approach to solve the GTSP exactly is 
a branch- and-cut algorithm [5(. By using this algorithm, 
Fischetti et al. solved several instances of size up to 89 
clusters; solving larger instances to optimality is still too 
hard nowadays. Two approximation algorithms for special 
cases of the GTSP were proposed in the literature; alas, the 
guaranteed solution quality of these algorithms is rather 
low for the real-world applications, see [2] and references 
therein. 

In order to obtain good (but not necessarily exact) so- 
lutions for larger GTSP instances, one should consider the 
heuristic approach. Several construction heuristics, dis- 
cussed in [2j, lZ|, [l4 1 , generally produce low quality solutions. 



A range of local searches, providing significant quality im- 
provement over the construction heuristics, are thoroughly 
discussed in 
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An ejection chain algorithm exploiting 
the idea of the TSP Lin-Kernighan heuristic is successfully 
applied to the GTSP in [10(. Although such complicated 
algorithms are able to approach the optimal solution by 
only several percent in less than a second for relatively 
large instances (the largest instance included in the test 
bed in [l(J has 1084 nodes and 217 clusters), higher quality 
solutions may be required in practice. In order to achieve a 
very high quality, one can use the metaheuristic approach. 
Among the most powerful heuristics for the GTSP, there is 
a number of memetic algorithms, see, e.g. , @, 0, S 
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Several other metaheuristic approaches were also applied 
to the GTSP in the literature, see, e.g., [II [lj [lj] . 

In this paper, we focus on a metaheuristic approach 
called ant colony optimization (ACO). ACO was first in- 
troduced by Dorigo et al. Q to solve discrete optimization 
problems and was inspired by the real ants behaviour. Ob- 
serve that, even without being able to see the landscape, 
ants are capable of finding the shortest paths between the 
food and the nest. This becomes possible due to a spe- 
cial substance called pheromone. Roughly saying, an ant 
tends to use a path with the highest pheromone concentra- 
tion. At the beginning, there are no pheromone trails, and 
each ant walks randomly until it finds food. Then it heads 
to the nest leaving a pheromone trail as it walks. This 
pheromone trail makes this path attractive to the other 
ants, and so they also reach the food and walk to the nest 
leaving more pheromone along the path. 

An ant does not necessarily follow the pheromone trail 
precisely. It may randomly select some slightly different 
path. Now assume that there are several paths between 
the food and the nest. The shorter is the path, the more 
frequent will be the walks of the ants using this path and, 
hence, the more pheromone it will get. Since pheromone 
evaporates with time, longer paths tend to get forgotten 



while shorter paths tend to become popular. Thus, in the 
end, most of the ants will use the shortest path. A more 
detailed description of the logic staying behind the ACO 
algorithms can be found in [3| and [4|. 

Since ants are capable of finding the shortest paths, it is 
natural to model their behaviour to solve such problems as 
the TSP or the GTSP. Several metaheuristics exploiting 
the idea of the ant colony, are proposed in the literature. 
In this study, we focus on the Ant Colony System (ACS) 
as it is described in 

There are two ACO implementations for the GTSP 
presented in the literature. The first one is an ACS heuris- 
tic by Pintea et al. It is an adaptation of the TSP 
ACS, and its performance is comparable with the most 
successful heuristics proposed by the time of its publica- 
tion. The second implementation by Yang et al. [18| is a 
hybrid ACS heuristic featured with a simple local search 
improvement procedure. 

We propose a new hybrid implementation of the ACO 
algorithm for the GTSP. The main framework of the meta- 
heuristic is a straightforward modification of the 'classical' 
TSP ACS implementation extended by an efficient local 
search procedure. We show that such a simple heuristic 
is capable of reaching near-optimal solution for the GTSP 
instances of moderate to large sizes in a very limited time. 

The paper is organized as follows. In Section 2, we 
briefly present the details of the ACS algorithm for the 
TSP. In Section 3, we propose several modifications needed 
to adapt the TSP algorithm for the GTSP. In Section 4, 
we describe the local search improvement algorithm used 
in the metaheuristic, and in Section 5, we report and anal- 
yse the results of our computational experiments. The 
outcomes of the research are summarized in Section 6. 

2. Basic ACS algorithm 

In this section, we briefly present the 'classical' ACS 
algorithm as described in Q. It is described for the TSP 
defined by a node set V of size n and distances d uv for every 
pair u 7^ v € V. If w(T) is the weight of a Hamiltonian 
cycle T (also called tour), the objective of the problem is 
to find T that minimizes w(T). 

A hybrid ACS algorithm is a metaheuristic repeat- 
edly constructing solutions, improving them with the lo- 
cal search procedure and updating the pheromone trails 
accordingly, see Algorithm [T] 

Let K be the set of ants. The typical number of ants 
I if I is 10. Let T k be an ordered set of nodes corresponding 
to the path of the ant k 6 K and T k be the ith node in T . 
Note that if \T k \ = n, the set T k can be considered as a 
tour. Let T^est be the best tour known so far. Initially, we 
set Tb os t Tnn; where Tnn is the tour obtained with the 
Nearest Neighbor TSP heuristic, see, e.g, [(| for description 
and discussion. 

At the initialization phase, the ants are randomly dis- 
tributed between the nodes: T k — {v}, where v £ V is 
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Algorithm 1 A high-level scheme of the hybrid ACS al- 
gorithm. 

Initialize pheromone trails, 
while termination condition is not met do 
Construct ants solutions. 
Apply local pheromone update. 
Improve the ants solutions with the local search 

heuristic. 
Save the best solution found so far. 
Apply global pheromone update, 
end while 



selected randomly for each k £ K. An initial amount 
To = w (Tx N ) °f pheromone is assigned t uv <— tq to each 
arc (u,v) £ E. This amount has to prevent the system 
from a quick convergence but also should not make the 
convergence too slow. 

On every iteration, each ant constructs a feasible TSP 
tour, which takes n — 1 steps. Let A kt C V be the set of 
nodes that the ant k 6 K can visit on the ith step, t = 
1, 2, . . . , n— 1. Since, in the TSP, an ant can visit any node 
that it did not visit before, A kt = V \ {T?,Tg, . . . ,T k }. 
Let r] uv be the so called visibility calculated as r\ uv = ^ — . 
Let a uv — r uv {ri uv )P , where (3 is an algorithm parameter, 
be the value defining how much attractive is the arc (u, v) 
for an ant. With the probability qo (that is an algorithm 
parameter selected in the range < qo < 1), the ant k, 
located in the node u = T k , selects the node v £ A kt 
that maximizes a uv . Otherwise it selects the node v £ A kt 
randomly, where the probability of choosing v is 



P 



kt 



y 



kt Q>h 



(2) 



On every step of an ant k £ K , & local pheromone 
update is performed as follows: 



(1 - £)T m + 



n ■ w(T NN ) 



(3) 



where £ is an algorithm parameter selected in the range 
< £ < 1 . This update reduces the probability of visiting 
the arc uv by the other ants, i.e., increases the chances of 
exploration of the other paths. 

After Ti — l steps, each T k for k € K can be considered 
as a feasible TSP tour. Run the local search improvement 
procedure for every T k and update the tour T k accord- 
ingly. The typical local search improvement procedure 
used for the TSP is fc-opt for k = 2 or k = 3. Now let 
k' = argmhifcg/f w(T k ) be the ant that performed best 
among K in this iteration. If w(Tk>) < w(J\, est ), update 
the best tour Tb es t found so far with TV. 

Finally, perform the global pheromone update. In global 
pheromone update, both evaporation and pheromone de- 
posit are applied only to the edges in the best tour Tbest 
found so far. Let p be an algorithm parameter called evap- 
oration rate and selected in the range < p < 1. Then 



the global pheromone update is applied as follows: 

P 



(1 - p)r uv + 



w{T hcst ) 



for (u, v) e T bost . (4) 



Before proceeding to the next iteration, reinitialize T k 
with {i>}, where v £ V is selected randomly for every k £ 
K. 

Various termination conditions can be used in an ACS 
algorithm. The most typical approaches are to limit the 
running time of the algorithm or to limit the number of 
consequent iterations in which no improvement to the orig- 
inal solution was found. 

3. Algorithm modifications 

In order to adapt the ACS algorithm for the GTSP, we 
need to introduce several changes. 

1. The Nearest Neighbor algorithm is redefined. Let T" 
for v £ V be a GTSP tour obtained as follows. Let A 
be a set of nodes. Set A <- V \ C{v). Set T v <- {v}. 
On every step t = 1, 2, . . . , m — 1, set T£ +l <— u and 
A <— A \ C(u), where u £ A is selected to minimize 
w(Tt,u). The output Tnn of the Nearest Neighbor 
heuristic is the shortest tour among T v , v £ V. 

2. The number \K \ of ants in the system is taken as an 
algorithm parameter and is discussed in Section [S] 

3. Since a GTSP tour visits only m nodes, the number 
of steps needed for an ant to construct a feasible tour 
is to — 1 . 

4. The set A kt of the nodes available for the ant k at 
the step t is defined as 



A kt = V \ \J C{T k ) . 



Let T^ est be the best tour found on or before the ith. 
iteration. The termination criteria used in our implemen- 
tation is as follows: terminate the algorithm if j > A and 
^ est = T hcst for i = j - A + l,j - A + 2, . . . ,j, where j 
is the index of the current iteration and A is an algorithm 
parameter. 



4. Local Search Improvement Heuristic 

It was noticed that many metaheuristics such as genetic 
algorithms or ant colony systems benefit from improving 
every candidate solution with a local search improvement 
procedure, see [l2| and references therein. Observe that 
all the successful GTSP metaheuristics are, in fact, hybrid. 
Thus, it is important to select an appropriate local search 
procedure in order to achieve a high performance. 

An extensive study of the GTSP local search algo- 
rithms can be found in 



11| . According to the classifi- 



cation provided there, all the local search neighborhoods 
considered in the literature can be split into three classes, 



3 



namely 'Cluster Optimization' (CO), 'TSP-inspired' and 
'Fragment Optimization'. While the latter one needs addi- 
tional research in order to be applied efficiently, neighbor- 
hoods of the other two classes are widely and successfully 
used in the metaheuristics, see, e.g., @, H, 15, 16| . 

The CO neighborhood is the most noticeable neighbor- 
hood in the CO class. Being of an exponential size, it can 
be explored in the polynomial time. Let T = (T±, T2, . . . , T m ) 
be the given tour. Then the CO neighborhood Nco(T) is 
defined as 

N C o(T) = {(TlX,---,Tj: 

T[ eC(Ti) for i = 1,2,. ..,m}. (5) 

Note that the size of the CO neighborhood is 

m 

\N co (T)\ = Y[\C l \£0(s m ), 

1=1 

where s = max™ 1 \Ci\ is the size of the largest cluster in 
the problem instance. Next we will briefly explain the CO 
algorithm finding the shortest tour T" £ Nco{T). 

Let T = (T U T 2 , . . . ,T m ) be the given tour. Create 
a copy S of the cluster C(Ti). Construct a multilayer di- 
rected graph G C o(Vco, E C o) with the layers C(T X ), G(T 2 ), 
. . . , C(T m ), S. For every pair of consecutive layers L\ and 
L2, for every pair of vertices ugli and v £ L2, create an 
arc (u,v) of weight d uv . Let P v be the shortest path from 
v £ C(Ti) to its copy v' £ S. Note that P v corresponds 
to a tour visiting the clusters in the same order as T does. 
Select v £ C(Ti) that minimizes the weight of P v . The 
corresponding cycle is the shortest tour T' £ Nqo (T) , and 
the procedure terminates in 0{ns 2 ) time. 

Several heuristic improvements of the above algorithm 
were proposed fill ] . In this research, we implemented only 
the easiest and the most important one. Note that the 
complexity of the algorithm linearly depends on the size 
of the cluster C(Ti). Since a tour can be arbitrarily ro- 
tated, let C(Ti) be the smallest cluster. This modifica- 
tion reduces the time complexity of the CO algorithm to 
0(njs), where 7 = min™]^ \C\\ is the size of the smallest 
cluster. 

Recall that the most typical neighborhoods used for the 
TSP are fc-opt. Several adaptation of the TSP fc-opt were 
proposed in [11], and the resulting neighborhoods were 
classified as 'TSP-inspired'. Since we aim at designing a 
fast and simple metaheuristic, we chose the 'Basic' fc-opt 
adaptation [ll|. In short, let Pgtsp be the original GTSP 
and let T = (Ti, T 2 , . . . , T m ) be the given tour defined in 
T^gtsp- Let Gtsp(VtsPi -Etsp) be the complete subgraph 
of G, where Vtsp = T 2 , . . . , T m }. Construct a TSP 
T'tsp for the graph Gtsp- Note that the tour T defined 
for Pqtsp is a feasible tour of the same weight in Vtsp, 
and any feasible tour in 'Ptsp is a feasible tour of the same 
weight in "Pgtsp ■ Improve the tour T with the TSP fc-opt 



algorithm. The obtained tour is the result of the 'Basic' 
adaptation of the fc-opt local search. 

It was shown that a combination of neighborhoods of 
different classes is often superior to the component local 
searches @. Thus, we use a local search that combines the 
neighborhoods of the CO and the 'TSP-inspired' classes. 
In particular, the improvement procedure used in our algo- 
rithm proceeds as follows. First, the given tour is improved 
with the 'Basic' adaptation of the 3-opt local search. Then, 
the CO algorithm is applied to it. No further optimiza- 
tion is performed so that the resulting solution is not guar- 
anteed to be a local minimum with regards to the 3-opt 
neighborhood. 

This local search procedure was obtained empirically 
after extensive computational experiments with different 
local search neighborhoods and strategies. 

5. Computational Experiments 

As a part of our research, we conducted extensive com- 
putational experiments to find the best parameter values 
and to measure the algorithm's performance. Our testbed 
includes a number of instances produced from the standard 
TSP benchmark instances by applying a simple clustering 
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We selected the 
Our ACS algo- 



by many researchers, see, e.£ 
same set of instances as in 
rithm and the local search procedures are implemented in 
C# and the computational platform is based on 2.93 GHz 
Intel Core 2 Due CPU. 

We used the following values of the algorithm param- 
eters: (3 = 3, p = 0.4, £ = 0.03, q = 0, A = 300 and 
\K I = 10. Among all the combinations of (3, p, £, <?o, A 
and \K\ that we tried, this one provided the best, on aver- 
age, experimental results. However, we noticed that slight 
variations of these values do not significantly change the 
behaviour of the metaheuristic. 

The extension of the local search procedure with the 
CO algorithm is the most significant modification imple- 
mented in our ACS. Thus, we start from studying the im- 
pact of the CO algorithm on the performance of the ACS. 
In our first series of experiments, we show the importance 
of this modification. In what follows, HACS refers to our 
hybrid ACS metaheuristic with the composite local search 
procedure as described above, and HACSo refers to the 
simplified version of the metaheuristic that uses only the 
3-opt algorithm as the local search procedure. 

The HACS and the HACSo algorithms are compared 
in Table [TJ The columns of the table are as follows: 

1. 'Instance' is the name of the the GTSP test instance. 
It consists of three parts, namely the number of clus- 
ters m, the type of the instance (derived from the 
original TSP instance) and the number of vertices n. 

2. 'Best' is the objective of the best solution known 
so far for the given problem instance. For the in- 
stances of size m < 89 the optimal solutions are 
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Table 1: Comparison of the HACS algorithm with its simplified version HACSo- 



known, see [5[. For the other instances the values 
are taken from Q. 

3. 'Error' is the relative solution error e, in percent, 
calculated as follows: 

e= »Cr)-«(3U) . 1Q0% 

w(T hcst ) 

where T is the solution to be evaluated and Tb cs t is 
the best solution known so far. 

4. 'Time' is the running time of the algorithm. 

5. 'Optimal' is the number of runs, in percent, in which 
the best known so far solution was obtained. 

The best result in a row is underlined. Since the ACO 
algorithms are non-deterministic, in order to get some sta- 
tistically significant results we repeat every experiment 10 
times. Hence, every result reported in Table [1] is an aver- 
age over the 10 runs. 

It is easy to see that the full version of the HACS 
clearly dominates the simplified one. This shows the im- 
portance of selecting the optimal nodes within clusters and 
also proves the efficiency of the approach used in our local 



search improvement procedure. It is worth noting that a 
more common adaptation of a TSP local search for the 
GTSP is to hybridize the 'TSP-inspired' and 'Cluster Op- 
timization' neighborhoods [ll, 14 1. However, our experi- 
ments prove that applying two local searches of different 
classes one after another may be a more effective strategy. 

In order to evaluate the efficiency of the HACS, we 
compare its performance to the performance of several 
other metaheuristics, see Table [5J In particular, we com- 
pare the HACS to three other metaheuristics, namely the 
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memetic algorithm SG by Silberholz and Golden 
memetic algorithm BAF by Bontoux et el. Q and an ACO 
algorithm PPC by Pintea et el. [ljjj]. 

The running times of SG and BAF reported in Tabled 
are normalized to compensate the difference in the exper- 
imental platforms. The SG algorithm was implemented in 
Java and tested on a machine with 3 GHz Intel Pentium 
4 CPU which we estimate to be approximately 1.5 times 
slower than our platform. The BAF algorithm was imple- 
mented in C++ and tested on a machine with 2 GHz Intel 
Pentium 4 CPU which we estimate to be similar to our 
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1.09 


10 


.81 


45ts225 


0.01 


0.14 


0.04 


0.03 


2.71 


1.14 


31 


.45 


46pr226 


0.00 


0.00 


0.00 


0.03 


2.29 


1.03 


8 


.25 


53gil262 


0.41 


0.45 


0.14 


0.22 


5.57 


2.43 


24 


.34 


53pr264 


0.00 


0.00 


0.00 


0.00 


3.83 


1.57 


18 


.27 


60pr299 


0.03 


0.05 


0.00 


0.24 


5.98 


3.06 


21 


.25 


641in318 


0.00 


0.00 


0.00 


0.12 


4.87 


5.39 


26 


.33 


80rd400 


0.62 


0.58 


0.42 


0.87 


9.95 


9.72 


32 


.21 


84A417 


0.00 


0.04 


0.00 


0.57 


7.22 


5.43 


31 


.63 


88pr439 


0.00 


0.00 


0.00 


0.78 


10.06 


12.71 


42 


.55 


89pcb442 


0.09 


0.01 


0.19 


0.69 


13.41 


15.62 


62 


.53 


99d493 


0.51 


0.47 


0.44 




22.68 


23.81 


166 


.10 


107att532 


0.15 


0.35 


0.05 




17.82 


21.13 


137 


.54 


107si535 


0.02 


0.08 


0.07 




19.99 


17.57 


90 


.98 


113pa561 


0.13 


1.50 


0.42 




19.26 


14.05 


1 49 


4S 


115rat575 


1.52 


1.12 


1.16 


— 


26.79 


32.32 


157 


.01 


131p654 


0.00 


0.29 


0.01 


— 


18.57 


21.78 


144 


.95 


132d657 


0.21 


0.45 


0.30 


— 


37.43 


88.16 


259 


.11 


1 A Kn'70A 
140U ( Z4 


1.0 ( 


U.O ( 


i no 




A Q QPI 
4o.oU 


1U f .GO 


218 


.66 


157rat783 


1.37 


1.17 


1.10 




47.41 


101.43 


391 


.79 


201prl002 


0.28 


0.24 


0.27 




123.38 


309.57 


513 


.48 


207sil032 


0.37 


0.37 


0.11 




177.00 


161.58 


616 


.28 


212ul060 


0.66 


2.25 


1.31 




103.89 


396.43 


762 


.86 


217vml084 


0.66 


0.90 


0.64 




95.35 


374.69 


583 


.44 


Average (all) 


0.33 


0.43 


0.30 




32.00 


66.61 


173 


.92 


Average (m < 89) 


0.09 


0.10 


0.06 


0.27 


5.66 


4.72 


25 


.40 



Table 2: Comparison of the HACS algorithm with the other GTSP metaheuristics. 



platform (note that the CH — h implementations are often 
considered to be twice faster than the Java or C# imple- 
mentations Q). The running time of PPC for each of the 
instances is 10 minutes as this was the termination criteria 
chosen in [l3[ (the computational platform is not reported 
in [H). 

For all the SG, BAF and PPC algorithms, the reported 
values are the averages among 5 runs; the results of HACS 
are the averages among 10 runs. 

Since the results of the PPC algorithm are reported for 
only a subset of the instances in our testbed, we provide 
two averages in every column of Tabled The first average 
(denoted as 'all') is the average over all the instances in our 
testbed, i.e., 40 < m < 217. The second average (denoted 
as m < 89) is the average over the testbed chosen in fl3j |. 
i.e., 40 < m < 89. 

In fact, we also compared our HACS to the ACO algo- 
rithm YSML by Yang et el. [18| and a memetic algorithm 
GK by Gutin and Karapetyan [7j , though those results are 
excluded from Table [5] 



The results reported in [18j are obtained for the in- 
stances of size 10 < m < 40 (the testbed was generated 
from the TSP instances by using the same clustering pro- 
cedure). It was noticed that these instances are relatively 
easy to solve to optimality even with a local search pro- 
cedure, see [l(j. Our ACS also solves all these instance 
to optimality and takes at most 1 sec for each run. The 
running time of YSML is not reported in [l8| , but the solu- 
tions obtained in [l8| are often not optimal. We conclude 
that our algorithm outperforms YSML. 

The GK memetic algorithm [7] is the state-of-the-art 
algorithm that, until now, was not outperformed by any 
other metahcuristic. It is a sophisticated heuristic with 
a well-tuned local search improvement procedure and in- 
novative genetic operators. Although GK dominates the 
HACS with respect to both the solution quality and the 
running time, it does not affect the outcomes of our re- 
search. Indeed, we aim at showing that a simple modifica- 
tion of the 'classical' ACO algorithm can yield an efficient 
solver for a hard combinatorial optimization problem. Also 
note that HACS and GK belong to the different classes of 



6 



metaheuristics. 

Table [2] shows that our HACS algorithm is similar to 
SG and BAF and significantly outperforms PPC with re- 
gards to the solution quality. Although, on average, BAF 
performs slightly better than HACS, there is no clear dom- 
ination since for some instances the HACS produces better 
solutions than BAF does. Similarly, SG is dominated by 
neither HACS nor BAF. With regards to the running time, 
HACS is the fastest heuristic for the large instances while 
SG usually takes less time for the instances of size m < 84. 
The BAF algorithm is the slowest one in every experiment 
and, on average, it is 5 times slower than HACS. 

Note that the above comparison of the running times 
is rather inaccurate since the considered algorithms were 
tested on different platforms, and only a rough normal- 
ization of the running times was performed. Still, certain 
outcomes can be made. In particular, the SG algorithm 
performs very well for the small instances while it is out- 
performed by HACS for larger instances with regards to 
both the solution quality and the running time. BAF, 
on average, produces better solutions then either HACS 
or SG do but this is achieved at the cost of significantly 
larger running times. Finally, HACS is superior to the 
other ACO algorithms, namely PPC and YSML, though 
the comparison was only possible for a limited number of 
test instances. 

6. Conclusions 

An efficient ACO heuristic for the GTSP is proposed 
in this paper. It is obtained from a 'classical' TSP ACS al- 
gorithm by several straightforward modifications and hy- 
bridisation with a simple local search procedure. It was 
shown that, among other reasons, the success of our HACS 
is due to the effective combination of two local search 
heuristics of different classes. Extensive computational 
experiments were conducted in order to prove that HACS 
performs as well as the most successful memetic algorithms 
proposed for the GTSP with the exception of the state- 
of-the-art sophisticated metaheuristic. It was also shown 
that HACS outperforms two other ACO GTSP algorithms 
proposed in the literature. 
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